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ABSTRACT 

We revisit the process of gravitational sedimentation of helium and heavy elements 
in the intra-cluster medium. We find that helium applies an inward drag force on 
heavy elements, boosting their sedimentation speed to nearly half its own. This speed 
is almost independent of the mass and the electric charge of heavy elements. In the 
absence of small-scale magnetic fields, helium sedimentation can increase the He/H 
abundance ratio in the cores of hot clusters by three orders of magnitude. It also 
steepens the baryonic density profile yielding a higher X-ray luminosity, which offers 
an explanation of the observed luminosity-temperature relation. 

If the primordial He/H ratio is assumed, then the gas density inferred from the 
observed X-ray emissivity might be underestimated by 30% in the cores of clusters 
and overestimated by 7% in the outer regions. The dark matter density on the other 
hand might be overestimated by a factor of 8/3 in the cores and underestimated by 
18% in outer regions. 
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1 INTRODUCTION 

In the equilibrium state of a multi-component plasma the 
number density, ni, of particles of mass mi, is n\ oc 



where 



is the gravitational potential. In 



the high temperature intra-cluster medium (ICM), light el- 
ements, helium (He) and hydrogen (H), can diffuse fast 
enough to reach their equilibrium distribution in a Hubble 
time. Larger frictional drag forces work on heavier elements, 
but still we expect, at least partial segregation as a result of 
diffusion. 

Diffusion can have important consequences. Fabian & 
Pringle (1977) suggested diffusion as a possible explana- 
tion for the observed gradients in the iron abundance inside 
galaxy clusters. In their calculations iron ions sediment by 
a distance comparable to the radius of the cluster within a 
Hubble time. For other elements they predicted the diffu- 
sion velocity to be proportional to AZ~ 2 , where A is the 
atomic number and eZ is the charge of the ion. These cal- 
culations assumed that only protons apply appreciable drag 
forces on heavy elements. Rephaeli (1978) claimed that by 
neglecting helium drag Fabian & Pringle overestimated the 
iron diffusion rate. 

Helium sedimentation can potentially change the global 
observational properties of X-ray emission from rich clusters. 
As has been noted by Qin & Wu (2000) current estimates 



of cluster masses from X-ray observations, which rely on the 
assumption of constant mean molecular weight, can be off 
by ~ 18% if helium is concentrated inside the core. Here 
we also point out that by steepening the baryonic density 
profile helium sedimentation increases the X-ray emissivity 
in hot clusters. 

In this paper we revisit the calculation of element sed- 
imentation in the ICM. We confirm the claim by Rephaeli 
that helium drag on iron is comparable to that of protons. 
However, instead of hindering the sedimentation of iron and 
other heavy elements we show that helium acts as a catalyst. 
Our analysis includes electric fields that are inevitably gen- 
erated by segregation of charged elements. Although these 
fields reduce the diffusion rate, the sedimentation time-scales 
of heavy elements are still several times shorter than previ- 
ous estimates. 

The paper is organized as follows. In §2 we present the 
basic equations and estimate the sedimentation speeds and 
time-scales. In §3 we discuss the equilibrium distribution as 
a limiting case of element sedimentation. We conclude in §4. 



© 0000 RAS 



2 Chuzhoy & Nusser 



2 THE EQUATIONS OF ELEMENT 
DIFFUSION 

We write the equations governing the evolution of individ- 
ual species in an ICM of any composition. We do not con- 
sider magnetic fields in this paper, but include electric fields 
which must exist in any ionized plasma in a gravitational 
field (Eddington f926). Let the ICM be made of any num- 
ber of species each made of particles characterized by mass 
mi and electric charge qi = eZ\. We denote the local num- 
ber density and velocity of a patch of matter of each species 
by m and Vi. Then, the mass density is p\ — nimi and the 
pressure is Pi = riikT, where we have assumed that the ICM 
is in local thermal equilibrium so all species share the same 
temperature T. With this notation the equations are 



dVj 

dt 



+ Vi ■ Wi 



VP qiE 
+9+ — 

pi rm 



where E is the electric field, and g is the gravitational field 
(force per unit mass). The constants Tlj are the time-scales 
for the drag forces acting on the species i as a result of 
Coulomb interactions with species j. Momentum conserva- 
tion implies that 7\j = (pj/pi) T ii- H mi 3> rrij then ry can be 
approximated by (Spitzer 1968) 



3(27r) 1/2 (fcT) 3/2 rm 
8irZ?Zfe 4 nj\nA m }/ 2 



(2) 



9.4 x 10' 
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where the InA is the Coulomb logarithm. The equations (1) 
must be supplied by additional relations that specify the 
electric field. To create an electric force comparable with 
gravitational and pressure forces a tiny fractional charge ex- 
cess (Gm p /e 2 ~ 10 -36 ) is sufficient. The additional relations 
can then be obtained by assuming charge neutrality and zero 
electric currents, i.e., 



^2 ni * = > 



and 



^2 "i<7i^i = . 



(3) 



2.1 Sedimentation speeds and time-scales 

Multiplying the equations (1) by n;m.i / ^ nimi and sum- 
ming over all species yields 

VP 

P 



9 



(4) 



where p — ^ nimi is the total local density of all species, 
V = "Y^jpiVi/p and P = ^2 Pi are, respectively, the mass- 
weighted mean velocity and total pressure, and Q = g — 
dV /dt — V ■ VV is the gravitational field in the frame of 
reference moving with the velocity V. The terms involving 
the electric and drag forces have disappeared by virtue of 
charge neutrality and momentum conservation, respectively. 
Since |Vi —V\ is much smaller than the thermal velocities, 
the velocity dispersion term ^Pi(Vi — V)V( V, — V)/p is 



negligible compared to the pressure and gravity terms we 
have not included it in the equation (4). 

At the initial stage of the diffusion process the dom- 
inant light species have the same distribution, which 
gives V(nikT)/(mmi) = (fj,m p /m,i)§, where /i = 
(l/m p ) Y2 norn/ ^2 rii is the mean molecular weight. Thus 
equations (1) yield 



( M m p /mi _ ^ + J^iVi - Vi)/ ni = 



(•>) 



Using equations (|3|) and (5) we find that the electric field is 
E = pg/e. From (2) 7ij cx nj' 1 rn^ 1 , so because of the low 
abundance of heavy ions and the low mass of electrons, only 
the drag terms from protons and helium ions is important. 
So we are left with a single equation for each element 



'1 + Zi 
.^47 



P 



1 9 + 



V p -Vi 







(6) 



We can now estimate the velocity of each species rel- 
ative to V p , the velocity of the proton fluid at each point. 
The advantage of using velocities relative to protons is the 
independence of the results from physical processes other 
than diffusion (radiative cooling, heating by supernovae etc). 
Even though these processes can have important dynamical 
effects on the cluster, they do not affect the expressions we 
develop for the relative velocities and element abundance. 



From (6) the helium velocity is 
V a -V p = (3/i/4 - l)g r ap w -0.56ff r ap . 



(7) 



The relative sedimentation speed of heavy elements, 
Vi — V p , can easily be related to that of helium V a — V p . 
Heavy elements experience an upward proton drag and in- 
ward helium drag forces. According to (2), Ty oc Z^ 2 and for 
large Z x these drag terms dominate all others in (6). Thus 
we are left with 

(V p -Vi)/T ip ^(Vi-V a )/n a . (8) 
Substituting n a and n p from (2) in the last equation yields 
V.-Vp 



V 



V 8n a J 



0.4, 



0) 



where we have taken n a /n p = 0.08. Thus, contrary to the 
prediction of Fabian & Pringle (1977), who by neglecting 
helium drag obtained (Vp — Vi) oc AiZ^~ 2 , the diffusion speeds 
are approximately the same for all heavy elements. Including 
all forces in the calculation increases slightly this result (by 
a factor of ~ 1 +4/Zi) 



V„ V, 



= 0. 



, 3Ai , 

1 + 



1 

~Z~i 



i 

~ 2 



(10) 



There is also a correction of a similar magnitude if the initial 
distribution of heavy elements is different from that of light 
elements. 

The relation (10) is independent of the physical state of 
the ICM. To obtain the relative velocities we have to know 
the temperature, the density and the gravitational acceler- 
ation. Taking in (2) InA ~ 40 as the typical value for the 
ICM, we find the helium velocity relative to protons to be 
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V p = 5 x lO 4 ^^- 1 T 8 



(11) 



where g 



5/(10" 



2 ), n p _ = n p /(10 3 m~ 3 ), and 



T s = T/10 8 K. This is lower by ~ 40% than the result 
of Qin & Wu (2000), the main difference is the inclusion 
of electric field in our calculation. As seen from equation 
( p"l| ) the diffusion speeds depend on the local density and 
temperature. However, a single diffusion time-scale is ob- 
tained if the gas and dark matter both follow an isothermal 
spherical density profile p(R) oc R~ 2 and the gas is in hy- 
drostatic equilibrium (g — g) with constant T. In this case 
Vi — Vj oc g/n p oc R, which together with continuity equa- 
tion, dm/dt = -R- 2 d R {R 2 Vi), yields, 



dt \ nt 



3(Vj ~ Vi) 
R 



(12) 



This motivates us to define the time-scale for diffusion be- 
tween two species as 



R 



m-v s y 



which for helium relative to protons gives 



where /b is the baryonic mass fraction in the cluster. 



3/2 



(13) 



(14) 



3 THE EQUILIBRIUM DISTRIBUTION 

The diffusion time-scale of helium as seen from (14) is com- 
parable to the Hubble time. It is thus prudent to examine in 
detail the equilibrium distribution of hydrogen and helium, 
which is the final product of diffusion. 

We assume a spherically symmetric cluster in which the 
dark mass inside a radius r is given by (Navarro, Frenk & 
White 1997) 



M dla (r) =4-K Ps rj 



ln(l + r/r s ) + 



1 + r/r s 



(15) 



where p s and r s are constants. The gas density profile, p, is 
determined by the equation of hydrostatic equilibrium (see 
eq.4), 



kT d^p 



/im p dr 



GM d Jr) 



(16) 



where we have neglected the contribution of baryons to grav- 
ity and assumed constant temperature, T, throughout the 
ICM. In the absence of sedimentation p is constant and the 
solution to (16) is (Makino, Sasati, & Suto 1998) 



p(r) = p(0)e-^(l + r/r s )^ r ° /r 



(17) 



where r\ — 4irGm p p s r 2 /kT. Sedimentation introduces a de- 
pendence of p on r and the above analytical solution is 
no longer valid. The abundances of the various elements 
are then determined by the equation of hydrostatic equi- 
librium for each element separately and the condition for 



local charge neutrality, 

kT d(lnrii) GM d 
nii dr r- 

ni(r)qi = . 



(18) 
(19) 



Taking m a = 4m p and neglecting the electron mass, these 
equations have the solution 



n p = 6Ci [/(r) + rV)-l] /h\r) 
n a = Ci[/(r)+/- 1 (r)-l] 2 /^M 



where f(r) — 



C 2 h 5 (r) - 1 + y/(C2h?(r) - l) 2 



(20) 
(21) 

1/3 



h(r) = (1 + r/r s ) r)rs,r , and C\ and C* 2 constants. The be- 
havior of the solution in the inner and outer regions is easily 
understood as follows. In the inner regions helium is domi- 
nant, thus p — 4/3 and E — pg / e — 4g/3e. From this follows 
rip tx (l+r/r s )~ v/s andn a oc (l+r/r s ) 4,|/3 . Note that, since 
the total force felt by protons becomes repulsive (eE > mg), 
their density is falling towards the center. The outer regions 
consist almost entirely of hydrogen plasma, thus p — 1/2 
and E = pg/e = g/2e. This gives n p oc (1 + r/r 3 ) v ^ 2 and 

oc rip, in agreement with Gil- 



n a oc (1 + r/r B ) iv y and so n a 
fanov & Syunyaev (1984). 



The constants C\ and C2 in the analytic solution (21) 
are fixed by the boundary conditions imposed on the abun- 
dances. Here we require that the ratio of total helium to 
hydrogen abundances inside the virial radius is equal to the 
primordial value (~ 0.08). We will present results for a clus- 
ter with gas temperature such that rj = 10 and a virial 
radius equal to 3r s in agreement with observations and N- 
body simulations of massive clusters (NFW 1997; Ettori & 
Fabian 1999). In Fig.[j] we show as the solid line the baryonic 
density obtained from the analytic solution. For comparison 
we also plot as the dashed curve the profile (17) which cor- 
responds to equilibrium without diffusion. A proper estima- 
tion of the density profile from the observed X-ray emissivity 
(oc (J2 niZ 2 ^2 n-iZi)) should take into account variations of 
the He/H abundance ratio throughout the ICM. Assuming 
constant He/H ratio can yield a biased estimate of the pro- 
file. To demonstrate this we plot as the dotted line in fig. [I] 
the estimated profile if a constant abundance ratio were as- 
sumed. In fig. |2| we show the baryonic mass fraction as a 
function of radius for the same three cases as in fig. jj]. We 
see (solid line) that diffusion introduces distinct features in 
the behavior of the baryonic fraction as a function of ra- 
dius. Finally fig. [] shows the number density of helium and 
hydrogen in the case with diffusion. 



4 CONCLUSIONS 

We have seen that diffusion steepens the gas density pro- 
file near the center of hot clusters, increasing their X-ray 
luminosity (by factor of 5 in our example). In colder clus- 
ters (T ~ 10 7 K) the diffusion time-scale is larger than the 
Hubble time and their luminosity remains unchanged. This 
may explain the observed discrepancy between the observed 
L oc T 3 (Mushotzky 1984; Edge & Stuwart 1991; David et 
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Figure 1. Gas density profiles. The solid and dashed lines corre- 
spond to equilibria with and without diffusion, respectively. The 
dotted line corresponds to the profile inferred from the X-ray 
emissivity by assuming constant abundances in the case with dif- 
fusion. The density is normalized so that the total baryonic mass 
within 3r a is unity. 




1 2 10 1 r/r 10° 



Figure 2. The same as the previous figure, but for the baryonic 
mass fraction (divided by the primordial value). 



al. 1993) and L oc T 2 which is expected from self-similar 
arguments (e.g., Kaiser 1986). 

If the inner regions of clusters are dominated by he- 
lium then the baryonic mass density as inferred from the 
X-ray emissivity can be underestimated by ~ 30% if con- 
stant He/H abundance ratio is assumed, while in the helium- 
deficient outer region it can be overestimated by ~ 7%. Esti- 
mates of of the dark matter density can be affected by even 
a larger factor. These estimates assume hydrostatic equilib- 
rium (eq. 16), so by taking fi — 0.59 (cosmic abundance) in- 
stead of 0.5 (pure hydrogen plasma) in the outer regions, we 
underestimate the total mass by ~ 18% (Qin & Wu 2000). In 
the helium dominated core the mass would be overestimated 
by a factor of 2.3. 

Our estimates of the sedimentation time-scales have 
neglected magnetic fields. Magnetic fields with coherence 
length comparable to the size of the cluster force the ions 
to move on longer orbits defined by the field lines. This 
can increase the sedimentation time-scales by factor of a 
few. Small-scale magnetic fields, however, can increase the 
time-scales by a factor ranging from a few in some estimates 
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Figure 3. The number density of hydrogen (solid) and helium 
(dashed) as a function of radius. The curves are normalized so 
that the total number of particles inside 3r s is unity. 

(Narayan & Medvedev 2001; Malyshkin 2001) to 100-1000 
in others (Chandran & Cowley 1998). 
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